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By  S .  G.  Robson 


ABSTRACT 

.Di  i sopropyl'methylphosphonate  (DIMP)  is  an  organic  compound  produced  as  a 
byproduct  of  manufacture  and  detoxification  of  a  GB  nerve  gas  agent  that  ap¬ 
pears  to  function  as  a  nonreactive  tracer  in  the  ground-water  environment. 

A  combination  of  ground-water  and  surface-water  transport  has  spread  tho  com¬ 
pound  through  large  areas  of  a  shallow  alluvial  aquifer  near  the  Rocky  Moun¬ 
tain  Arsenal. 

A  ground-water  solute  transport  model  has  been  used  to  simulate  DIMP  con¬ 
centrations  in  the  alluvial  aquifer.  The  Phase  I  model  is  the  first  and  small¬ 
er  of  two  models  to  be  built  in  order  to  evaluate  the  effectiveness  of  various 
ground-water  quality  management  alternatives.  The  model  covers  an  area  of 
about  30  square  miles  (77.7  square  kilometres)  and  was  calibrated  primarily 
by  comparison  of  simulated  and  observed  1975  DIMP  concentrations.  Two  prelim¬ 
inary  20-year  model  projection  runs  were  made  to  demonstrate  the  predictive 
ability  of  the  model.  The  results  of  the  larger  Phase  II  model  and  the  final 
model  projections  will  be  presented  in  a  forthcoming  final  report. 


INTRODUCTION 


Historical  Background  of  the  Contamination  Problem 

The  Rocky  Mountain  Arsenal  was  established  in  19^2  and  presently  encom¬ 
passes  about  27  mi^  (70  kni‘)  immediately  northeast  of  the  city  of  Denver,  Colo. 
The  northern  boundary  of  the  arsenal  is  about  7  mi  (11  km)  south  of  the  city 
of  Brighton,  Colo.  The  arsenal  was  constructed  to  produce  toxic  chemical  and 
Incendiary  munitions — about  2k2,000  tons  (220,000  t)  of  which  v/ere  manufac¬ 
tured  during  World  War  II  (Reynolds,  1975).  After  a  postwar  reduction  in  ac¬ 
tivities,  the  arsenal  began  in  1953  a  second  manufacturing  phase,  upon  comple¬ 
tion  of  a  new  toxic  chemical  agent  facility.  This  manufacturing  and  filling 
facility  manufactured  a  GB  nerve  gas  agent  from  1953  to  1957-  One  byproduct 
of  this  manufacturing  process  was  the  organic  compound  d i i sopropy Imethy 1 phos- 
phonate  (DIMP).  The  arsenal  has  subsequently  been  involved  in  manufacturing 
anticrop  agents  and  removing  toxic  chemicals  from  munitions.  From  1973  to  the 
present  (1975)  the  GB  manufacturing  facility  has  been  used  to  detoxify  the  GB 
agent  held  in  bulk  storage  and  in  bombs  and  warheads. 
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The  manufacturing  and  detoxification  activities  at  the  arsenal  have  pro¬ 
duced  industrial  effluents,  some  of  which  have  high  di ssol ved-sol ids  concen¬ 
tration  and  (or)  high  concentrations  of  organic  compounds.  Prior  to  1956 
these  effluents  were  disposed  of  in  unlined  surface  ponds.  Effluent  percolat¬ 
ing  through  the  bottom  of  these  ponds  has  degraded  the  chemical  quality  of  the 
ground  water  in  a  shallow  alluvial  aquifer.  Subsequent  to  1956  industrial  ef¬ 
fluents  have  been  discharged  to  a  lined  reservoir  (Reservoir  F)  and  the  unlin¬ 
ed  ponds  have  been  either  unused  or  used  temporarily  to  hold  surface  runoff  or 
fresh  water. 

Petri  and  Smith  (1956)  studied  the  ground-water  quality  near  the  arsenal 
in  1955  and  1956.  They  found  that  high  concentrations  of  sodium  and  chloride 
in  the  aquifer  extended  from  near  the  arsenal  disposal  ponds  northwest  to  the 
South  Platte  River.  Data  from  this  and  subsequent  work  were  used  by  Konikow 
(1975a)  to  prepare  detailed  maps  of  the  bedrock  surface  elevation,  v/ater-table 
elevation,  and  aquifer  saturated  thickness  and  transmissivity  near  the  arsenal. 
Konikow  later  expanded  this  work  to  include  a  digital  ground-water-quality 
model  capable  of  simulating  the  movement  and  dispersion  of  chloride  in  the 
aquifer  (Konikow,  l?75b). 

In  197^1  OIMP  was  detected  in  ground  water  on  the  arsenal  and  surround¬ 
ing  property  to  the  north  and  west.  Trace  concentrations  v/ere  detected  in 
wells  as  much  as  7  mi  (11  km)  downgradient  from  the  arsenal  disposal  ponds  and 
within  1.0  mi  (1.6  km)  of  two  wells  in  the  city  of  Brighton's  municipal  well 
field  (pi.  l) .  Concern  has  buen  expressed  about  OIMP  in  potable  ground  vrater 
because  little  is  known  about  toxicity  levels  or  the  effect  of  long-term  con¬ 
sumption  of  trice  amounts  of  the  compound.  This  and  an  associated  problem  of 
ground-water  contamination  by  a  pesticide  byproduct  led  the  Colorado  Depart¬ 
ment  of  Health  to  issue  Cease  and  Desist  Orders  against  the  Rocky  Mountain 
Arsenal  and  the  Shell  Chemical  Co.  The  Shell  Chemical  Co.  plant  on  the  ar¬ 
senal  property  is  though'  to  be  the  source  of  the  ground-v/atcr  contamination 
by  the  pesticide  byproduct  d icyclopentad icne  (OCPD).  The  Cease  and  Desist 
Orders  state,  in  general,  that  the  Rocky  Mountain  Arsenal  and  the  Shell  Chem¬ 
ical  Co.  are  to: 

1.  Immediately  stop  the  off-post  discharge  of  DIMP  and  DCPO  from  both 
surface-  and  ground-water  flcxi. 

2.  Submit  a  proposed  plan  of  action  to  preclude  such  future  off-post 
discharge  and  take  action  on  the  approved  plan. 

3.  Develop  and  institute  a  surveillance  plan  to  verify  compile  ice  with 
items  1  and  2. 


Purpose  and  Scope 


The  Rocky  Mountain  Arsenal  requested  that  the  U.S,  Gcolog i ca 1'  Survey  con¬ 
duct  a  d  i  g i ta 1 -mode  1  study  of  the  DIMP  contamination  of  ground  'water  on  the 
arsenal  and  surrounding  property  in  order  to  provide  insight  into  the  mechan- 
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Ism  of  DIMP  transport  and  to  provide  data  useful  in  meeting  the  requirements 
of  Cease  and  Desist  Orders.  In  particular,  information  is  needed  concerning 
the  area  affected,  concentrations,  and  rate  and  direction  of  movement  of  the 
plume  of  contaminated  ground  v/ater.  If  Brighton's  municipal  well  field  (pl.l) 
will  be  affected,  the  model  may  be  used  to  determine  future  arrival  times  and 
concentrations. 

The  study  is  being  undertaken  in  two  phases.  Phase  I  involves  the  con¬ 
version  of  the  existing  chloride  transport  model  (Konikow,  1975b)  of  the  ar¬ 
senal  to  a  OIMP  transport  model  in  order  to  evaluate  the  magnitude  of  the 
ground-water  contamination  in  the  immediate  area  of  the  arsenal.  Phase  II 
involves  building  a  second  DIMP  transport  model  of  a  larger  area  In  order  to 
investigate  the  potential  for  DIMP  moving  into  Brighton's  well  field  and  other 
adjacent  areas.  Both  models  will  assume  conservative  (nonreactive)  transport 
of  DIMP  and  wiil  be  calibrated  using  present  data  on  DIMP  concentrations  in 
the  aquifer.  This  report  documents  the  results  of  Phase  I.  A  subsequent  re¬ 
port  will  document  the  work  done  on  Phase  M. 


Acknowledgments 

Existing  data  on  the  OIMP  concentrations  in  ground  and  surface  waters 
were  furnished  by  the  Rocky  Mountain  Arsenal.  These  data  included  analyses 
run  by  the  Colorado  Department  of  Health,  Shell  Chemical  Co.,  and  the  Rocky 
Mountain  Arsenal.  Tlie  arsenal  aiso  ran  O'MP  analyses  on  supplemental  samples 
collected  by  the  U.S.  Geological  Survey. 


MODEL  DEVELOPMENT 

The  model  used  for  both  the  chloride  and  DIMP  transport  studies  at  the 
arsenal  is  based  on  an  iterative  alternating-direction  implicit  mathematical 
solution  of  the  ground-water  flow  equations  coupled  with  a  method  of  charac¬ 
teristics  solution  of  the  solute  transport  equations  as  described  by  Brede- 
hoeft  and  Pinder  (1973).  These  mathematicai  procedures  require  that  the  area 
modeled  be  divided  into  numerous  small  rectangular  segments  or  nodes  of  equal 
dimensions.  In  the  Phase  I  model  a  1 ,000-ft  (300-m)  interval  grid  was  used 
to  delineate  the  nodes  (pi.  l).  At  each  node  in  the  area  to  be  modeled  the 
geohydrologic  and  chemical  characteristics  of  the  aquife'"  are  described. 

The  geohydrologic  characteristics  of  the  alluvial  aquifer  near  the  arsen¬ 
al  have  been  described  by  Konikow  (1975a).  In  gei.eral,  he  found  that  the  di¬ 
rection  of  ground-water  movement  is  from  the  southeast  to  northwest  tov/ard  the 
South  Platte  River  (pi.  l),  and  that  the  saturated  thickness  of  the  alluvial 
aquifer  ranges  from  0  to  60  ft  (0  to  18  m) .  The  aquifer  is  discontinuous  2  to 
3  mi  (3  to  5  km)  east  of  the  South  Platte  River  because  numerous  bedrock  highs 
separate  areas  in  which  the  alluvium  is  saturated.  In  this  area  the  direction 
of  ground-water  movement  is  controlled  by  the  configuration  of  the  aqui/er. 

The  transmissivity  of  the  aquifer  near  Reservoir  F  is  about  3,000  ft,^/d  (200 
m^/d),  v/hich  is  about  7  times  lower  than  the  transmissivity  near  the'  river 
(fig.  1). 
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Water-level  hydrographs  for  wells  north  and  west  of  the  arsenal  show  that 
there  have  been  minimal  long-term  water-level  changes  in  the  area  betv/een  1937 
and  1975  (Brookman,  1969) •  The  most  significant  short-term  water-level  change 
occurred  during  the  195^”57  drought  when  water  levels  declined  3  to  5  ft  (0-9 
to  1.5  m)  .  For  modeling  purposes,  the  c'-ound-wotcr  system  was  assumed  to  be 
in  a  steady-flow  condition  from  the  time  'f  DIMP  discharge  (1952-56)  to  the 
present  (1975).  To  simulate  these  conditions,  a  steady-state  flov/  model  was 
coupled  with  a  transient-state  transport  model.  This  modeling  procedure  ena¬ 
bled  the  simulation  of  the  steady  ground-water  flow  system  in  conjunction  v;ith 
aguifer  DIMP  concentrations  that  change  during  the  1952-75  simulation  period. 

The  boundary  conditions  for  the  DIMP  model  were  handled  in  the  same  man¬ 
ner  as  were  those  for  the  chloride  model  (Konikow,  1975b).  Boundaries  trav¬ 
ersing  the  aquifer  were  modeled  as  constant-head  boundaries;  that  is,  bound¬ 
aries  at  which  the  rate  of  flow  across  the  boundary  may  vary  but  the  head  at 
the  boundary  is  held  constant.  No-flov/  boundaries  were  simulated  between  the 
aquifer  and  the  bedrock  outcrops.  Although  the  bedrock  is  known  to  contain 
permeable  zones  (McCcnaghy  and  others,  1964),  the  hydraulic  conductivity  of 
these  zones  is  generally  much  less  than  that  of  the  alluvium.  As  a  result, 
the  assumption  of  impermeable  model  boundaries  is  thought  to  be  valid.  This 
assumption  is  further  substantiated  by  the  successful  model  simulation. 

The  quantity  and  distribution  of  recharge  and  discharge  to  the  DiMP  model 
were  similar  to  those  used  In  the  chloride  model  (pi.  1  and  table  1),  with 
one  exception.  Industrial  effluent  high  in  chloride  was  disposed  into  ponds  A 
through  E  from  about  1943  to  1956.  Effluent  containing  DIMP  was  thought  to  be 
disposed  Into  ponds  A  through  E  from  1952  to  1956.  Unfortunately,  no  informa¬ 
tion  is  available  on  the  concentration  of  DIMP  in  the  effluent.  As  a  result, 
the  DIMP  source  concentrations  used  in  the  model  were  determined  by  trial  and 
error  as  those  which  produced  the  best  agreement  between  the  1975  field  data 
and  the  corresponding  model  calculations. 

The  model  assumes  that  conservative  (nonreactive)  ground-water  transport 
and  dispersion  of  DIMP  v/i 1 1  occur  as  the  contaminant  moves  through  the  aqui¬ 
fer.  In  the  absence  of  test  data,  it  can  only  be  assumed  that  DIMP  is  conser¬ 
vative;  that  is,  it  does  not  adsorb  on  the  soil  clay  matrix,  combine  with  oth¬ 
er  compounds,  or  vindergo  oxidation  or  reduction  as  it  moves  through  the  aqui¬ 
fer.  This  assumption  appears  to  be  valid  v^hen  the  concentration  distributions 
of  DiMP  and  the  conservative  tracer  chloride  are  compared  and  is  further  sub¬ 
stantiated  by  the  fact  that  the  same  transverse  and  longitudinal  dispersivlty 
(100  ft  or  30  m)  proved  adequate  in  both  the  chloride  and  DIMP  models. 


MODEL  CALIBRATION 

A  calibration  procedure  was  used  to  check  the  validity  of  the  DiMP  trans¬ 
port  model.  Calibration  involved  a  comparison  of  (l)  The  observed  water-table 
elevations  with  model -calculated  water-table  elevations,  and  (2)  observed  mid- 
1975  DIMP  concentrations  with  the  model -calculated  1975  DIMP  concentrations. 
The  agreement  between  the  observed  and  calculated  quantities  is  indicative  of 
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Table 

1 . — Recharge, 

discharge , 
for  Phase  I 

arid  concentre 
DU'fP  model 

tion  of  DIt-iP  data 

Source' 

1952-75  average 

recharge  or  discharge, 

in  cubic  feet  per  second  in  m 

Average  01  HP 

concentrat i on , 

i crograms  per  litre 

Recharge 

Discharge 

1952-56 

1956-7l(  197I1-75 

A 

0.072 

- 

17,000 

0  0 

B 

.012 

- 

17,000 

0  0 

C 

.662 

- 

17,000 

0  3,000 

D 

.059 

- 

1 ,000 

0  0 

E 

.027 

- 

1  ,000 

0  0 

F 

0 

0 

0 

0  c 

1 

6.528 

- 

* 

*  A 

L 

.01(9 

- 

0 

0  0 

N 

.075 

- 

0 

0  0 

0 

2.162 

- 

0 

*  * 

R 

- 

11.813 

- 

- 

S 

- 

.  1.258 

- 

- 

u 

1(.667 

- 

0 

0  0 

w 

- 

1 .21(2 

- 

- 

'See  plate  1  for  location  of  sources. 

*0lMP  concentration  in  recharge  varies  during  1952"75  period. 
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the  ability  of  the  model  to  simulate  either  historical  or  future  water-table 
elevations  and  OIMP  concentrations  in  the  aquifer.  The  model-calculated 
water-table  elevations  are  generally  within  3  ft  (0.6  m)  of  those  shown  on 
plate  1,  and  all  the  model -generated  water-level  contours  are  located  within 
700  ft  (210  m)  of  the  corresponding  contour  shown  on  plate  1.  The  similarity 
between  the  observed  DIMP  concentrations  and  the  calculated  concentrations  may 
be  seen  by  comparing  figures  2  and  3*  The  lack  of  historical  data  on  the  quan¬ 
tities  of  water  discharged  to  each  pond,  the  concentration  of  DIMP  in  this  wa¬ 
ter,  and  the  concentrations  of  DIMP  in  the  ground  water  prior  to  1375  contri¬ 
bute  uncertainties  to  the  model  which  are  evidenced  as  discrepancies  between 
the  calculated  and  the  observed  data.  The  agreement  attained  between  the  mod¬ 
el  calculations  and  the  observed  data  is  sufficient,  however,  to  consider  the 
model  calibrated  and  valid. 

Because  ground-water  OIMP  concentration  data  arc  not  available  for  years 
prior  to  197^,  the  model  calibration  is  based  on  known  concentrations  at  the 
beginning  and  end  of  the  simulation  period.  The  lack  of  data  for  Intermediate 
checks  on  calibration  is  not  a  severe  limitation  because  of  the  congruence  be¬ 
tween  the  DIMP  transport  model  and  the  chloride  transport  model.  Chloride 
data  were  available  for  intermediate  calibration  checks  on  the  chloride  model. 
As  the  hydrologic  factors  which  affect  chloride  transport  also  affect  DIMP 
transport,  the  proper  calibration  of  the  chloride  model  adds  support  to  the 
less  extensive  calibration  of  the  DIMP  model. 

If  further  OiMP-model  calibration  is  desired,  adequate  concentration 
changes  probably  will  have  occurred  in  the  aquifer  after  2  to  3  years.  A 
short-term  calibration  would  then  be  possible  if  adequate  data  are  collected 
in  the  interim.  These  data  need  to  include  (1)  DIMP  concentrations  in  the 
sources  of  ground-water  recharge  and  in  the  aquifer,  (2)  the  quantity  of 
ground-v;ater  recharge  from  ponds  and  other  sources,  (3)  water-level  measure¬ 
ments  in  wells,  and  (A)  the  quantity  of  water  pumped  from  wells  and.  sumps. 

The  calibrated  DIMP  model  was  used  to  calculate  the  ground-water  flow 
velocities  in  the  alluvial  aquifer  (fig.  ^) .  Velocities  ranged  from  less  than 
0.1  ft/d  (0.03  m/d)  in  the  southeast  part  of  the  model  area  to  more  than  I5 
ft/d  (0.A6  m/d)  near  the  South  Platte  River.  When  figure  k  is  used  in  con¬ 
junction  with  the  direction  of  ground-water  movement  shown  on  plate  1,  the 
ground-water  travel  time  between  adjacent  points  in  the  aquifer  may  be  cal¬ 
culated. 

The  model  calculates  annual  DIMP  concentrations  in  the  ground  water  re¬ 
sulting  from  the  estimated  average  concentration  and  rate  of  recharge  from 
each  pond  from  1952  to  1375.  By  examining  the  spatial  distribution  of  DIMP  in 
the  aquifer  at  different  times  it  is  possible  to  better  understand  the  mechan¬ 
ism  by  which  the  DIMP  contamination  reached  its  present  locations. 

The  DIMP  concentrations  in  the  aquifer  were  assumed  to  be  zero  in  1952, 
prior  to  the  disposal  of  DIMP  to  the  environment.  During  the  following  L 
years,  DIMP  was  disposed  of  in  the  unllned  ponds  A,  B,  C,  D,  and  E '(pi ,  l) . 
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Figure  3.— Model -generated  1975  concentration  of  DIMP  in  ground  water. 


OF  EQUAL  GROUND-WATER  VELOCITY, 
FEET  PER  DAY— Literval  Is  variable, 
convert  feet  squared  per  day  to  metres 
larcd  per  d.ay  multiply  by  0.30^3 


Figure  k. — Model-calculated  ground-water  velocities 
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The  resulting  model  calculations  show  DIMP  concentrations  in  excess  of  1,000 
yg/1  (micrograms  per  litre)  in  the  areas  shown  on  figure  5-  After  the  com¬ 
pletion  of  the  lined  Reservoir  F  in  1956,  it  is  assumed  that  ground-v/ater  re¬ 
charge  of  D I MP-contaminated  water  ceased.  The  resulting  model  calculations 
show  an  improvement  in  ground-vjater  quality  near  ponds  C,  D,  and  E  (fig.  5). 
The  tv^o  zones  in  excess  of  1 ,000  ug/1  near  ponds  A  and  B  had  merged  by  I960 
and  remained  near  this  location  through  subsequent  years  due  to  the  low 
ground-water  velocities  in  the  area.  By  examining  figures  5  through  9  it  can 
be  seen  that  the  main  part  of  the  contaminated  ground  water  moves  from  ponds 
C,  D,  and  E  to  the  South  Platte  River  through  the  aquifer  to  the  v/est  of  Res¬ 
ervoir  F.  The  model  indicates  that  by  1972  (fig.  9)  this  plume  was  no  longer 
present  in  concentrations  greater  than  1,000  ug/l.  However,  DIMP  concentra¬ 
tions  in  excess  of  1,000  yg/1  persist  from  1956  to  1975  northeast  of  Reservoir 
F  due  to  the  lower  ground-water  velocities  in  this  area. 

Ground-water  discharge  occurs  at  springs  and  seeps  along  a  reach  of  First 
Creek  from  near  96th  Avenue  to  the  O'Brian  Canal.  The  model  indicated  that  a 
maximum  DIMP  concentration  of  about  2  yg/1  would  occur  in  the  canal  as  a  re¬ 
sult  of  First  Creek  water  discharging  into  the  canal.  This  is  consistent  with 
the  2.5  pg/I  DIMP  found  in  two  samples  of  water  from  Barr  Lake  (the  terminus 
of  the  O' Brian  Canal). 

When  no  additional  DIMP  input  after  1956  is  assumed  in  the  model,  the 
1975  calculations  show  an  area  larger  than  1  mi^  (3  km^)  near  ponds  C,  D,  and 
E,  and  Reservoir  F  with  a  DIMP  concentration  less  than  5  ug/1.  The  DIMP  con¬ 
centrations  observed  in  May  and  June  1975  in  this  area  range  between  about  300 
and  3>000  yg/1.  This  discrepancy  suggests  that  an  additional  input  of  DIMP 
has  occurred  in  this  area.  The  change  in  DIMP  concentration  with  time,  shov/n 
in  figure  10,  indicates  that  ground  water  with  a  high  DIMP  concentration  has 
moved  Into  the  area  of  v/ell  1*1  during  the  period  of  record,  August  197^  to 
August  1975*  (Locations  of  wells  discussed  are  shov/n  on  fig.  2.)  The  high- 
concentration  water  moved  into  well  1i*1  prior  to  the  start  of  the  record  in 
January  1975.  This  could  indicate  that  well  lAl  is  closer  to  the  source  of 
contamination  than  is  v/ell  41.  The  concentrations  found  in  v/ell  127  suggest 
that:  (l)  The  contaminant  arrived  at  this  well  at  about  the  same  time  as  at 

well  4l ,  and  (2)  v/ell  127  is  located  near  the  lateral  edge  of  the  degraded 
zone.  The  relatively  low  concentrations  of  DIMP  found  in  wells  145  and  129 
(fig.  2)  indicate  that  ground  v/ater  v/ith  DIMP  concentrations  as  high  as  those 
in  wells  4l  and  l4l  does  not  appear  to  be  moving  downgradient  from  ponds  A  or 
B.  The  direction  of  ground-water  movement  shov/n  on  plate  1,  in  addition  to 
the  above  data,  suggest  that  the  source  of  the  recent  DIMP  input  is  in  the 
southern  one-half  of  sec.  26,  T.  2  S.,  R.  67  W.  There  are  insufficient  data 
to  determine  the  mechanism  responsible  for  the  possible  recent  DIMP  input  or 
the  exact  location  of  the  source.  Possibilities  to  be  considered  include: 

1.  Percolation  in  ponds  C,  D,  or  E  of  surface  runoff  from  ponds  A  or  B. 

2.  Induced  recharge  of  Dl MP-contaminated  water  from  permeable  strata  in 

the  underlying  Dawson  Formation.  .  * 
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3.  Release  of  OIMP  previously  concentrated  by  evaporation  from  the  soil 
of  the  unsaturated  zone. 

A  leak  in  the  liner  or  discharge  pipeline  to  Reservoir  F;  however, 
there  are  no  conclusive  water-quality  data  to  indicate  that  a  leak  in  the  lin¬ 
er  of  Reservoir  F  exists..  It  was  found  by  use  of  the  model  that  a  DIMP  input 
from  pond  C,  of  1-year  duration  at  a  concentration  of  3.000  ug/1,  v/as  adequate 
to  reasonably  simulate  the  1975  observed  concentrations  in  the  aquifer  near 
ponds  C,  D,  and  E,  and  Reservoir  F. 


PRELIMINARV  MODEL  PROJECTIONS 

The  model  Is  able,  to  calculate  future  DIMP  concentrations  that  would  oc¬ 
cur  in  the  aquifer  as  a  result  of  a  given  set  of  assumptions.  The  model  was 
used  to  calculate  the  DIMP  concentrations  that  would  occur  during  the  next  20 
years  if  no  significant  changes  occurred  in  the  aquifer  system  or  the  opera¬ 
tional  practices  of  the  ponds.  Model  assumptions  v/ere  that:  (1)  The  future 
quantities  of  water  recharged  from  ponds  A  through  E  would  not  differ  from  the 
historic  average,  (2)  the  OIMP  concentration  of  this  water  would  be  zero,  (3) 
the  aquifer  would  remain  in  a  steady-flow  condition  In  the  future,  and  (L)  the 
pattern  and  quantity  of  irrigation  recharge  and  pumpage  would  not  change  over 
the  next-  20  years.  The  observed  1975  OIMP  concentrations  in  the  aquifer  (fig, 
2)  were  used  as  the  Initial  concentration  of  the  model  aquifer  for  1975.  The 
resulting  mode  1 -ca 1 cu 1 ated  concentrations  for  1985  and  1995  are  shown  on  fig¬ 
ures  11  and  12.  By  comparing  figures  2,  11,  and  12  the  changes  in  OIMP  con¬ 
centration  between  1975.  1985,  and  1995  may  be  seen.  The  model  indicates  that 
the  high  OIMP  concentrations  near  pond  C  in  1975  will  have  an  adverse  effect 
on  the  concentrations  in  the  aquifer  between  pond  C  and  the  South  Platte  River 
In  1985,  OIMP  concentrations  in  excess  of  20  ug/1  probably  will  occur  in  this 
area.  By  1995  this  zone  of  increased  concentration  probably  will  have  decreas 
ed  in  size  and  no  longer  extend  to  the  South  Platte  River,  A  general  trend  of 
improving  v/ater  quality  is  evident  in  most  of  the  remaining  area  of  contamina¬ 
tion  from  1975  to  1995. 

A  second  set  of  conditions  to  be  evaluated  by  use  of  the  model  involved 
the  installation  of  a  barrier  to  the  movement  of  ground  v/ater  between  ponds  B 
and  C,  The  model  assumptions  for  this  condition  were  the  same  as  before  ex¬ 
cept  that  an  impermeable  barrier  was  assumed  to  halt  all  underflo.-/  moving  past 
the  location  of  the  barrier.  No  water-level  changes  occurred  upgradient  of 
the  barrier  because  it  was  assumed  that  any  water  collecting  behind  the  bar¬ 
rier  would  be  artificially  removed  from  the  hydrologic  system.  Water  levels 
in  the  aquifer  immediately  downgradient  from  the  barrier  decli.ned  to  the  ex¬ 
tent  that  the  aquifer  was  dewatered.  The  resulting  10-  and  20-year  model  pro¬ 
jections  are  shown  on  figures  13  and  l1».  A  comparison  of  figures  2,  13,  and 
shows  the  changes  in  DIMP  concentrations  from  1975  to  1995.  The  high  con¬ 
centrations  near  pend  C  in  1975  still  have  an  adverse  effect  on  the  1985  con¬ 
centrations  in  the  aquifer  between  the  pond  and  the  South  Platte  Rrver.  The 
slower  improvenient  in  quality  in  this  area  by  1995  than  was  shov^n  on  figure  12 
is  due  to  the  reduced  ground-water  velocity  produced  by  the  loss  of  underflov/ 
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LINE  OF  EQUAL  DIMP  CONCENTRATION 
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past  the  barrier.  It  should  be  noied  that  the  reduction  of  underflow  through 
an  area  will  retard  the  movement  of  ground  water  but  will  do  little  to  change 
the  DIMP  concentration  of  the  water, 

MANAGEMENT  OBJECTVES 

Before  a  "best"  water-quality  management  procedure  can  be  determined,  the 
desired  objective  needs  to  be  firmly  established  by  the  Rocky  Mountain  Arsenal, 
the  Colorado  Department  of  Health,  and  other  concerned  parties.  This  objec¬ 
tive  will  strongly  influence  the  choice  of  management  procedures  to  be  used. 

Two  hypothetical  examples  can  be  used  to  illustrate  this  point.  In  the  first 
example,  the  objective  would  be  to  reduce  the  DIMP  concentrations  in  the  allu¬ 
vial  aquifer  both  on  and  off  the  arsenal  as  rapidly  as  possible.  In  the  sec¬ 
ond  example,  the  objective  would  be  to  reduce  the  amount  of  DIMP  reaching  the 
South  Platte  River  and  the  surrounding  alluvial  aquifer  outside  the  arsenal. 

A  water-quality  management  procedure  that  might  be  effective  in  meeting  the 
first  example  objective  would  be  to  Increase  the  quantity  of  fresh-water  re¬ 
charge  to  the  aquifer  on  the  arsenal.  This  (-ould  increase  the  g round -v/a ter 
velocities  in  the  aquifer  and  move  the  DlMP-contaminated  water  out  of  the  area 
more  rapidly.  However,  this  management  procedure  could  be  detrimental  if  the 
second  example  objective  were  to  be  used.  A  management  procedure  that  would 
be  better  suited  to  the  second  objective  would  be  to  cease  all  artificial  re¬ 
charge  on  the  arsenal  and  pump  numerous  wells  in  the  degraded  plume  near  the 
arsenal  boundary.  The  two  hypothetical  examples  requite  drastically  different 
management  procedures,  but  both  have  objectives  which  would  lead  to  an  im¬ 
provement  in  ground-water  quality. 

SUMMARY 

The  results  presented  in  this  report  complete  the  work  on  the  Phase  I 
DIMP  model.  The  model  has  been  developed  and  adequately  calibrated  by  com¬ 
paring  model  calculations  with  observed  field  data.  The  model  is  capable  of 
calculating  DIMP  concentrations  in  ground  v/ater  that  v/ould  occur  in  response 
to  a  variety  of  g-ound-woter  management  alternatives.  Two  examples  have  been 
given  to  illustrate  how  different  water-quality  management  procedures  can  af¬ 
fect  the  DIMP  concentrations  in  the  aquifer  over  a  20-year  period.  The  ex¬ 
amples  indicated  that  the  addition  of  a  barrier  to  the  movement  of  ground  wa¬ 
ter  located-  south  of  pond  C  would  reduce  the  quantity  of  water  moving  through 
the  downgradient  parts  of  the  aquifer,  thereby  reducing  the  ground-v^ater  ve¬ 
locities  in  these  areas.  As  a  result,  the  DIMP  concentrations  in  ground  water 
v/oiild  be  slightly  higher  at  the  end  of  20  years  with  the  barrier  than  they 
would  be  without  the  barrier.  Model  runs  using  different  management  altern¬ 
atives,  different  assumptions  about  the  ground-water  flow  sys.um,  or  different 
time  periods  are  also  possible.  The  exact  objective  of  the  ground-water  qual¬ 
ity  management  must  be  firmly  established  before  an  optimum  plan  of  action  can 
be  determined,  as  very  similar  objectives  can  require  drast'cally  different 
plans  of  action. 
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